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Abstract We consider the effect of parametric uncertainty on properties of Linear Time Invariant 
systems. Traditional approaches to this problem determine the worst-case gains of the system 
over the uncertainty set. Whilst such approaches are computationally tractable, the upper bound 
obtained is not necessarily informative in terms of assessing the influence of the parameters on 
the system performance. We present theoretical results that lead to simple, convex algorithms 
producing parametric bounds on the /i 2 -induced input-to-output and state-to-output gains as 
a function of the uncertain parameters. These bounds provide quantitative information about 
how the uncertainty affects the system. 
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1. INTRODUCTION 

The application of polynomial optimisation techniques to 
robust control problems have been successful both in terms 
of sum-of-squares programs (Blekherman et al. (2013)) 
and the generalized problem of moments (Lasserre (2009)). 
Some of these advances were highlighted in the survey 
paper on Robust Control (Petersen and Tempo, 2014, 
Section 6.1). The stability analysis of linear systems with 
linear dependence on uncertain parameters lying on the 
simplex was performed using parameter-dependent Lya¬ 
punov functions (LFs) as detailed in Chesi (2010) (see 
Section VI.C). In Chesi et al. (2005) it is proven that, 
for time-invariant systems, LFs which are quadratic in the 
state and homogeneous on parameters are necessary and 
sufficient to prove stability of the linear system. However 
the resulting test in terms of homogeneous polynomials 
may require a LF of large degree. Polynomial LFs have 
also been studied by Bliman (2004) where, exploiting the 
fact that the parameterised Lyapunov matrix is analytic 
on the uncertainty domain, it is demonstrated that there 
exists a homogenous polynomial LF proving the stability 
for uncertainties in the hypercube. In Oliveira and Peres 
(2007), Polya’s theorem is applied to solve homogeneous 
inequalities on the simplex. Conditions for discrete-time 
linear systems with polynomial parameter-dependent LFs 
were presented in Lavaei and Aghdam (2008) where the 


* Work supported by the Engineering and Physical Sciences Re¬ 
search Council projects EP/J012041/1. Giorgio Valmorbida is also a 
Fulford Junior Research fellow at Somerville College. James Ander¬ 
son acknowledges funding from St John’s College, Oxford. Dhruva 
Raman is supported by the EPSRC funded Systems Biology DTC 
at the University of Oxford. 


simplex as the uncertain domain is presented as a partic¬ 
ular case. 

Dissipation inequalities for uncertain linear systems were 
studied considering interval matrix uncertainties in Barb 
et al. (2003), where the limits on the uncertainty for 
which a dissipation inequality holds were investigated. 
In Sato and Peaucelle (2006), parameter-dependent LF 
were considered to compute bounds on Hoo and 'H 2 gains 
over polytopic uncertainty domains for linear systems with 
rational dependence on the uncertain parameters. The 
authors then use the slack variable method to solve the 
resulting set of matrix inequalities which are affine on the 
uncertain parameters. The same conditions for worst-case 
computation of "Hoc and 7^2 norms were also presented in 
(Chesi, 2010, Section VI.D). In the case of autonomous 
systems, £00 norm bounds on nonlinear system trajecto¬ 
ries, valid over compact sets of initial conditions, have been 
constructed using LF based approaches in Chesi and Hung 
(2008). 

A parameterised property (or parameterised bounds on a 
property), such as the 'Hoo-norm of an uncertain system 
gives the designer important information for the selection 
of the components of a physical system (which define 
the system parameters) and allows them to quantify be¬ 
forehand the degradation, measured as the norm of the 
mismatch to the target (nominal) performance. Predefined 
bounds on performance degradation can be useful for 
robust and resilient control law design. In this context, 
obtaining a closed-loop performance curve that respects 
a specified degradation bound becomes a requirement 
for control synthesis. This is in contrast to the custom¬ 
ary approach of searching for guaranteed cost solutions 
which accounts for the worst case scenario given a set 





of uncertain parameters. Also, in the context of gain¬ 
scheduling (Packard (1994)), having a parameterised per¬ 
formance index for a set of given controllers may be of use 
when designing the scheduling parameter. 

A system is said to be robust to variations in its parameters 
if it is able to maintain its performance properties subject 
to a bounded variation of these parameters. For instance, 
the difference between the output signal of the nominal 
model and that of any perturbed model is expected to 
be small in some sense. From the viewpoint of system 
identification, such a property can be problematic: a low 
degree of uncertainty in the input lends itself to a high 
degree of uncertainty in the parameter estimate, and the 
identification protocol becomes ill-conditioned. This prop¬ 
erty has been variously referred to as quantitative uniden- 
tifiability (Vajda et al. (1989)), practical unidentifiability 
(Holmberg and Ranta (1982)), and sloppiness (Brown and 
Sethna (2003)). 

Traditional methods of diagnosis (Vajda et al. (1989)) 
involve linearisation of the effects of parameter pertur¬ 
bation around a previously provided parameter estimate. 
However such effects are typically highly nonlinear (Hines 
et al. (2014)). Recent advances from various angles have 
been made in the construction of nonlinear approximations 
to regions of practically unidentifiable parameters around 
a nominal estimate (e.g. Raue et al. (2009); Calderhead 
and Girolami (2011)). Such advances typically rely on edu¬ 
cated sampling of outputs over parameter space, with each 
sample involving the simulation of an ODE. Parameterised 
performance indices as provided in this paper, in con¬ 
trast, provide algebraic characterisations of performance 
degradation in parameter space. This allows for simple 
construction of level sets of parameters inducing a given 
magnitude output perturbation, or optimisation of some 
system property subject to parametric constraints speci¬ 
fying a given maximal level of performance degradation. 


with a;(0) = xo, A : H : R"'’ R’^x'", 

C : R"'’ —>• RPX"- as the nominal system. Consider the 
linear time-invariant parameter-dependent system 

( X = A{9)x + B{0)u 
\y = C{e)i + D{e)u, 

a;(0) = a;(0) = Xq (we consider the initial conditions of (1) 
and (2) to match for reasons that will become clear in next 
section). System (2) is called the uncertain system. 


Define the signals e{t) := x{t) — x{t) and Ay(t) := y[t) — 
y{t) to obtain 
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( 3 ) 


where AA(6») := A{e*) - A(6»), AH(6») := B{9*) - B{9) 
and AC(6I) := C{9*) - C{9), AD{9) := D{9*) - D{9), 
a;(0) = Xo and e(0) = 0. 


We assume throughout that A(9) is Hurwitz V9 € 0. From 
the linearity of (3) and the triangular structure of the 
system matrix we have stability of the origin. 


In the following section we describe input and state-to- 
output properties of (3) considering the output Ay, which 
is the deviation from the nominal behavior. 


The output of the uncertain system (2) is written in terms 
of (x, e) as 


y=[C(9) -C(9)j 


D{9)u. 


The interest in introducing y as above, instead of the 
mismatched output Ay as in (3), is motivated by the fact 
that we may want a bound for ||?/(t )||2 — ||?/(i )||2 instead 
of \\y{t) — y{t )\\2 = ||Aj/(t )||2 as developed in this paper. 


This paper is organized as follows: Section 2 formulates 
the problem. In Section 3 we present results allowing to 
characterise input/state-to-output properties as functions 
of uncertain parameters. Section 4 illustrates the results 
with an analytical and numerical examples. The paper is 
concluded in Section 5. 

Notation. We take R" to denote the n-dimensional Eu¬ 
clidean space. R>o specifies the set of non-negative real 
numbers. / denotes the identity matrix. For a given matrix 
M S R"^", Tr{M) is the trace of M, det{M) is the deter¬ 
minant, and He{M) = M + M’^. For a symmetric matrix, 
if M > (<)0, then M is positive (negative) semidefinite. 
For Xq G R” and x{t) G £ 2 , ||a;o||, ||a :||2 respectively 
denote Euclidean norm of Xg and the squared £ 2 -norm 
of x(t). is defined to be the set of continuous, fc-times 
differentiable, scalar functions. 

2. PROBLEM FORMULATION 

Let 0 C R"*’ define the set of uncertain parameters. Let 
0* G 0 be the fixed nominal parameter values. We refer to 
the linear system 

( 1 ) 


2.1 Invariance of the output 

In order to make a distinction between the transient effects 
and the steady-state effects of the uncertainty on the 
mismatch signal Ay{t) let us define the set of parameters 
that leave the steady-state output invariant with constant 
inputs. 

Definition 1. (Steady-State Output Invariant Sets). The set 
in the space of parameters 

0« := (6» G R”" I lim Ay{t) = 0Vu, G r) (4) 

L t—^OO J 

is called the steady-state output invariant set of parameters 

with respect to i. Here, Ui G R"* is a constant input vector 

taking the value 0 in all components except the The 

set 

m 

:= n ©ol (5) 

i=l 

is then defined as the steady-state output invariant set of 
parameters. 

Let us now define the set of parameters leaving the output 
invariant with time-varying inputs 

Definition 2. (Output invariant set of parameters). The set 
0« := {9 G R"» I Ay{t) = 0Vt, u,{t)} (6) 


X = A{9*)x + B{9*)u 
y = Cl9*)x + D{9*)u 














is called the invariant set of parameters with respect to 
input i. Here, Ui is as defined in Definition 1. 

m 

©0 := n (7) 

i=l 

is the invariant set of parameters. 

Proposition 1. For system (3) the steady-state output 
invariant set of parameters as in Definition 1 with respect 
to inpnt i can equivalently be characterised by 

0 « :={ 0 eR"'> |/«.W= 0 } ( 8 ) 

with 

fss{e) = AC{e*)A-\e*)B,{0*) 

+ c{9)A-\e) [AB,{e) - AA-\e)A-\e*)B,{e*)]. 

Proof. Note that, for any constant scalar input up. 
i = 0 ^ Ay = 0, Ay = 0, V 6 l € 0 

Expand both sides of the implication using (3). Substitut¬ 
ing the left hand side equation into the two right hand side 
equations and simplifying gives the result. ■ 

By defining x = 

^ ( X = A{9)x + B{0)u 

^ ■ \y = C{e)x +d19)u. 

Proposition 2. For system (3) the output invariant set of 
parameters as in Definition 2 with respect to input i can 
equivalently be characterised by 

0^ := {9 € R"" I f,s{s,9) = OVs € C} (10) 

with 

fss{s, 9) = C{e*) {si - A{9)y^ y{9) + DA9). 

Proof. A parameter is in the output invariant set if and 
only if the following holds: 

Ay = 0 Ay = 0. 

This follows from the observation that Ay = 0 at time 
zero. Expanding both sides of the implication using (3) 
gives the result. ■ 

3. PARAMETERISED INPUT/STATE-OUTPUT 
PROPERTIES 

We propose a set of parameterised input/state-output 
properties for linear systems with the goal of character¬ 
izing the deviation from the nominal behavior. We assume 
that both nominal and uncertain system matrices are 
Hurwitz. No assumptions on the class of functions de¬ 
scribing the parameter dependence, or the uncertainty set, 
are required. A convex optimisation formulation to solve 
some dissipation inequalities introduced in this section is 
presented in Section 4, which restricts the dependence on 
the parameters to be of rational form. 

Proposition 3. (Upper Bounds on State-to-Output Gain). 
For any initial condition xq = x(0) and cq = 0, if there 
exists a function V : R" x R” x R”® —>• R>o, V € 

U( 0 , 0 ,-) = 0 and functions pi : R”® —>■ R>o, i = 1,2 and 
a set 0 C R”® such that 

V{x,e,9) <pi{9)\\x\\'^+P2{9)\\e\\'^ V 6 » € 0 (11a) 


y = 


V 

Ay 


we denote system (3) as 


V{x, e, 9) < -{Ayf{Ay) V 6 » £ 0 (11b) 

along the trajectories of (3) (with u = 0) then 

^^<Pi(0) V0e0. (lie) 

Proof. Integrate (11b) over the interval [0,T] to obtain 

f {Ayiyy{Ay{T))dT < V{x{0),e{0),e)-V{x{T),e{T),9). 
Jo 

Since V{x{T)^e{T),6) > 0 we have 

T 

(Ay(T))'^(Ay(T))dT < U(x(0), e(0), 6 >) 
for every T. With (11a) we obtain 

[ (Ay(r))'^(Ay(r))dr < pi{9)\\x\\^ + p2{9)\\e\\^. 

Jo 

Since the initial conditions of the system satisfy (x(0), e(0)) = 
(xo,0), and the relation holds for every T > 0 we ob¬ 
tain ( 11 c). ■ 

The above proposition helps in characterising sets of 
parameters bounding the £ 2 -norm of the mismatch output 
in terms of level sets of the function pi{9) since 

9g{9g 0|pi(0) <\}^ [ {Ayf{Ay)dT < A||xof. 

Jo 

In Section 4 an example illustrates level sets of a function 
Pi{9) satisfying the above conditions. 

Remark 1. Notice that we can always impose pi{9*) = 0 
when computing pi, p 2 and V to solve the inequalities 
(11). This fact results from the system structure and the 
choice of output to be the mismatch to the nominal system 
which implies that, whenever 9 = 9*, we have Ay{t) = 0 

(yielding || Ay||| = 0) which implies = 0. 

The proposition below provides conditions to lower bound 
the State-to-Output Gain of system (3). 

Proposition f. (Lower Bounds on State-to-Output Gain). 
For any initial condition xq = x(0), if there exists a 
function 1^ : R" X R" x R”® —)► R>o, Ve G and functions 
P£i : R”® —>■ R>o, * = 1,2 and a set 0 C R"® such that 

pn{9)\\xf +Pi2{9)\\e\\^ <Viix,e,9) V0 G 0 (12a) 

-{AyfiAy) < Vi{x, e, 9) V0 G 0 (12b) 

along the trajectories of (3) then 

^^>m( 0 ) V 0 G 0 . ( 12 c) 

Ikolr 


The following proposition proposes a parameterisation for 
the £ 2 -induced gain of system (3). 

Proposition 5. (Upper Bounds on the £ 2 -Induced Gain). 
For x(0) = 0, e(0) = 0, if there exists a function U : R" x 
RTi ^ R>o, U(0, 0,-) = 0, U G and functions 

: R"® —>• R>o and 7 „ : R"® — R>o and a set 0 C R"® 
such that 


V{x,e,9)>0 V6IG0 (13a) 

V{x,e,9) < - 7 d( 6 ')(Ay)^(Ay) +^n{9)u^u V 6 » G 0 

(13b) 

then 

llAyll? 


l«ll 


< 7(0) V0 G 0, 


2 

2 


(13c) 










with 7 ( 0 ) 


ln{0) 
ld{0) ' 


Proof. Integrate (13b) over the interval [0,T] to obtain 

J - 7 ( 6 I)(m(t))^m(t)) dr 

<F(x(O),e(O),0)-F(x(T),e(r),0). 
By hypothesis P(0, 0, •) = 0, which gives 

J {ld{d){Ay{T) f {Ay{T)) - 'yn{0){u{T))’^ u{t)) dr 

<-Vix{T),e{T),0). 

From (13a) we have 

/ {ld{d){Ay{T)f {Ay{T)) - 7 „( 6 »)(u(r))'^u(r)) dr < 0. 
Jo 


which is equivalent to (13c) since (13a)-(13b) hold for all 

6> e 0. ■ 

Remark 2. Whenever 0 is a compact set, an upper bound 
on the guaranteed T-Loc cost of the uncertain system is given 
by 7 *+maxe 7 (^), where 7 * is the "Hca gain of the nominal 
system. As pointed out in Remark 1, when solving (13), it 
is possible to impose 7 „( 6 >*) = 0 since for 6 = 9* we have 
Ay = 0 if a;(0) = e(0) = 0. 

The next result applies to system (9) with D{9) = 0. Recall 
that the 7 ^ 2 -norm of system ( 1 ) is 

\\G{9*)\\n, = TY[B^{e*)PB{d*)] = [C{9)QC^ (9*)] 

where P > 0 solves {9*)P + PA{9) + C'^{9*)C{9) = 0 
and Q > 0 solves the dual equation. An upper-bound 
for ||G( 0 *)|||^^ can be found by replacing the equality in 
the Lyapunov equation with an inequality and searching 
for any positive definite matrix P that solves the matrix 
inequality. 

Proposition 6. (Bounds on the P 2 -norm). If there exists a 
function R : R" X R" x R’"'’ R>o, R(0, 0, •) = 0, R e 
and functions qi,q 2 '■ Q ^ R>o satisfying, 

V[x,e,9) < qi{9)\\x\\^ + q2{9)\\e\\‘^ V 6 » e 0 (14a) 

V{x,e,9) + {Ay)'^Ay <Q Vd G 0 (14b) 

then 

\\G{9)\\n.<ymTr{B^{e)B{9)) V0 G 0. (14c) 

Proof. Integrating (14b) over the interval [0,T] and sub¬ 
stituting Ay{t) = C{9)x(t), C{9) := [AC{9)C{9)] gives 


V{x{T),e{T), 9) - fo(x(0), e(0), d) < - C {Ay{ t)^A y{r)dr 

^0 

= ~ll^y|l 2 ,[ 0 .T]- 

Letting t —>■ 00 we arrive at 

ll2/lll<mO),e(O),0). _ (15) 

Let G{s) be a frequency domain representation of G (note 
that for clarity we do drop the dependence on 9), then the 
'P 2 -norm is equivalently given by 


Define Y{t) := Ce^'^B and let yt denote the i**' column of 
Y, thus Pi = Ge^^Bi- Let ^^(t) = e^^Bi, then it follows 
that yi(r) = Cxi{T) and therefore 

m 

The bound (15) gives 

|| 2 /zf < fo(a;.(0), eo, 9) = R(e^"P„ e(0), 9) 

= R(P„e(O),0) 
<qi{9)BjB, (16) 

where the final inequality comes from e(0) = 0 and (14a). 
Finally, using (16) it follows that 

771 m 

E II 2 /*II < 91 ( 0 ) E = 9i(^)Tr[P^P]. 

which upon inserting the dependence on 9 completes the 
proof. ■ 

Corollary 7. (Bounds on the P 2 -norm). If there exists a 
parameter dependent matrix function P ; 0 —>■ R 2 nx 2 n 
and a function gi : 0 R>o satisfying, for all 0 G 0, 

0 < P{9) <qi{9)I (17a) 

A^{9)P{9) + P{9)A{9) + C{9fC{9) < 0 (17b) 

then 

\\G{9)f^^<qA9)Tr{B^{9)m G 0. (17c) 

4. EXAMPLES 
4-1 Analytic Example 
Consider system 

r. _ 9t 

E" 1 + 02 *'' (18) 

[y = x 

with 9* G R>o, i G {1,2}. The above system corresponds 
to (1) with A{9*) = -jSt, B{9*) = 0 and G{9*) = 1. 
Consider the uncertainty set 

0 = {d G R^|0i G (0,+oo)} 

since 0i > 0 and 1 + 02 > 0, V0 G 0, A{9) is Hurwitz for 
all 0 G 0. Now define 

0 '+=(^-i) *e{i, 2 } 

establishing a bijection between 0 and the set 
0 = G R^|0i G (-l,+oo)| . 


We can then construct system (3) for (24) in terms of 0 as 
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A{9*) 
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0{(l + 0i) 
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1 + 02(1 + 02) 
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l|G(s)||?,,= 


27r 


Tr [G*{ju:)G{ 3 oj)\ dio 


{Ce^^Bf{Ce^^B)dT 


In order to obtain bounds on state-to-output gain as a 
function of the parameters, following Proposition 3, let us 
consider the structure 


V{x,e,9) =p-^{9)x^ +P 2 ( 0 )e^ 


= Tr 


(19) 





















with Pi : Q ^ K>o- With the above, and using the fact 
that e = Ay, inequality (11b) is written as V (x, e, 9) + 

X 




< 0 with 


M(0) = 


/ 

Pii9)A{9*) p2{9)A{9*)^^f^^ y y 

\ 


1 + 9*{1 + 92) 



^ -P2i9)9t{l + 9i) 1 1 


V 

L 1 + 9*{1 + 92) 2 J 
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Let us now compute > 0 and P 2 > 0 to satisfy M{9) < 0 
V0 G 0 therefore verifying V (x, e, 0) + < 0 for all 0 G 0. 

First consider the element (2,2) in M(9*,9) 


^{2,2){9) — 


The choice of p 2 as 


2p2{e)ei{i + h) 

tT02(td”^2) 


P2{9) 


(1 + ^ 2(1 + ^ 2 )) 

9l{l+~9i) 


( 20 ) 


satisfies P2i9) > 0 V0 G 0 and yields M(^2,2)0) = ~1- With 
the above choice for p 2 we have 


detiMm = -p^mi9*)- 

( 21 ) 

replacing the value for P2{9) and picking 


Pi(^) 


1 91 nei{e2-h)-~ei) \ 

2(1 + 0 *)!^ 91 ( 1 + h) j 


( 22 ) 


we obtain det(M(9)) = 0. Since det(M{9)) = 0 and 
M(^ 2 , 2 ) < 0, the matrix M{9) in the quadratic form 
V (x, e, 9) + is negative semidefinite and we have 
V(x,e,9) + < 0V9 G 0. Notice pi(9) is non-negative 

V0 G 0 and Pi (9) = 0 in the subspace 


On 



Following Proposition 3 we have that 



Figure 1. Upper bound for 

in (22) for 0G [-1,0.5] x [-1,1]. 


given by pi(9) as 
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Figure 2. Level sets of pi(9). The set 0„ where pi(9) = 0 
is depicted by the dashed black line. The level set 
corresponding to the first conical region neighboring 
the dashed line is {9 j pi(9) < 1}. 


Pln(9) 


iMmi 

l|x(0)112 

(l|x(0)p) ^ /o°°(e^^® - e^(®)'^xo)^dr 

P g2A(e*)r _ 2 g(A(dn+A(e))T g2A(e)r^^ 

° 1 2 1 

~2A(9*) (A(9*) + A(9)) ~ 2A(9) 

(i+eU _ 2{i+ei){i+e^{i+02)) 

2^1 ( 0 * (1+^2 ( 1+^2 

i+e^ji+e^) 

2fl*(i+ei) • 


(23) 

The above expression therefore characterises the state-to- 
output gain of system (24) as a rational function of the 
parameters 9. Figure 3 depicts that pi in (22) is an upper 
bound to pi„. 


{Ay)'^(Ay)dT < V(x(0),0,9) =pi(0)x^(O), 

from the fact that pi(9) = 0 V0 G On we conclude that 
Ioi^yVi^y)dx = 0 V0 G 0„. 

We depict pi(9), the upper bound to IjAj/jjl, in Figure 1 
and level sets of pi{9) in Figure 2. 

For this simple linear scalar system it is actually possible 
to compute the £ 2 -norm of the output, and to obtain the 
state-to-output gain as: 



4-2 Numerical example 

In this example we make use of polynomial optimisation 
techniques to numerically construct parameter-dependent 
bounds on the £ 2 -induced norm of an uncertain system. 

Consider system 

J X = (—3 -\- 9i + 92 )x + (—1 -|- 02)u 
{y=(2-9i)x ^24) 

with 9* = (0, 0) and the uncertainty set defined as the rect¬ 
angle 0 = [—1.5,1.5] X [—1, Ij. In order to compute the C 2 - 
induced gain by checking the condition on Proposition 5 we 
















































2.6 



Figure 3. Surface Pin{S) as in (23) and its upper 
bound pi{0) as in ( 22 ). 

consider the following constraints (see (Blekherman et al., 
2013, p. 115) for further details) 

V{x,e,0) - mii{x,e,9)pei - mi2{x,e,9)pe^ G Y.[x,e,9] 

(25a) 

- V{x,e,9) - 7d(6>)(A?/)'^(Ay) +"fn{9)u^u 

- m 2 i{x, e, 9)pg^ - m 22 {x, e, 9)pe^ G E[a;, e, 9] (25b) 

mij{x,e,9) &Tj[x,e,9] i,_7'G{l,2} (25c) 

with polynomials = —(0i + l-5)(0i —1.5),p6i2 = ~(^i + 

l)(di — 1 ), which yields 0 = {d | > 0,^02 > 0 } to 

formulate the following Sum-of-Squares program 

minimize 7 „( 0 ) subject to (25). (26) 

We set polynomials V and to be quadratic on variables 
X and e and respectively of degree 3 and 2 on variable 9. We 
formulate (26) with SOSTOOLS (Papachristodoulou et al. 
(2013)) and solve the underlying SDP with SDPT3 (Toh 
et al. (1999)). Figure 4 depicts solutions to the convex 
optimisation Problem (26). The guaranteed cost 7 = 
over the set of uncertainties, is obtained by imposing 
deg(^n) = degipfd) = 0 in the dissipation inequality (25b) 
and is depicted by the light gray hyperplane. The darker 
gray hyperplane, which sits underneath the guaranteed 
cost, corresponds to an affine 7 „ {deg(^n) = 1 ) and a 
constant 7 ^. The gain curve that corresponds to the case 
deg(pin) = ‘^,deg{'yd) = 1 is a better approximation for 
the actual gain, depicted in the bottom graph, its level 
sets are depicted on the plane ( 01 , 02 )- Notice that there 
is no guarantee that the computed surfaces have empty 
intersection since they were computed independently. 

5. DISCUSSION AND CONCLUSIONS 

We presented conditions to obtain upper-bounds to 
input/state-to-output gains for uncertain linear systems. 
These upper-bounds are functions of the uncertain param¬ 
eters. For systems with polynomial and rational depen¬ 
dence on the parameters it is possible to compute these 
upper-bounds by solving a set of polynomial inequalities. 
These inequalities can then be solved with semidefinite 
programming. We also provide a single-dimensional exam¬ 
ple where the exact characterization of the state-to-output 
gain is possible and we compare this gain with an analytic 
solution to the proposed conditions. 


7(0) 



Figure 4. Graphs of the computed 7 ( 0 ) = for 

polynomials 7 ^( 0 ) and 7 d( 0 ) of different degrees. 
The depicted surfaces correspond to GC: Parameter- 
Independent Upper Bound, AC: Affine Parameter- 
Dependent Upper Bound, RC: Rational Parameter- 
Dependent Upper Bound, TC: Actual Cost. 

An extension of the formulation presented is the state- 
feedback control design with a given parameterised upper 
bound of the input-output gains. Such a parameterised 
upper bound represents an allowable degradation of the 
nominal trajectories. Another challenge is to assess the 
effects of uncertainty in the control gains while computing 
the control laws. Criteria for non-fragile control laws Keel 
and Bhattacharyya (1997) can be defined in terms of 
the upper-bounds on the system gains parameterised in 
terms of uncertainties in the control law gains. Quadratic 
separation and ellipsoidal sets for controller uncertainties 
were used in Peaucelle and Arzelier (2005) in this context. 

The extension to nonlinear polynomial systems (systems 
defined by polynomial vector fields) are currently under 
study. For this class further assumptions are required since 
the existence of a single equilibrium point may not be 
guaranteed with the parameter variation. One important 
aspect is that the cascaded structure of the system when 
described in terms of the original state x and the error e 
as in (3), is preserved. However, the fact that the stability 
of the origin of the elements of the cascade implies the 
stability of the origin of the overall system is no longer 
true if the system under analysis is non-linear. 
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Appendix A. SUM-OF-SQUARES POLYNOMIALS 

Let us take ]R[a;] to be the ring of polynomials with 
real coefficients in the vector x € H.” of indeterminates. 
Define the subring S[a;] be the set of sums of squares of 
polynomials in ]R[a:]: 

m 

E[a;] = G K[a;] | 3m G N, G K[a;] : p = 

2=1 


